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We investigate the phenomenological viability of a recently proposed class of composite dark mat¬ 
ter models where the relic density is determined by 3 —> 2 number-changing processes in the dark 
sector. Here the pions of the strongly interacting field theory constitute the dark matter particles. By 
performing a consistent next-to-leading and next-to-next-to-leading order chiral perturbative inves¬ 
tigation we demonstrate that the leading order analysis cannot be used to draw conclusions about 
the viability of the model. We further show that higher order corrections substantially increase the 
tension with phenomenological constraints challenging the viability of the simplest realisation of the 
strongly interacting massive particle (SIMP) paradigm. 
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I. INTRODUCTION 

Dark Matter accounts for circa 85% of fhe maffer in 
the universe, but besides from ifs cosmological abun¬ 
dance, very litfle is known abouf ifs nafure. In a 
wide class of models, fhe relic abundance is generafed 
via a thermal freeze-out in the early universe. T 5 rpi- 
cally 2 —> 2 annihilation processes into e.g. standard 
model particles keep dark matter in thermal equilib¬ 
rium with the standard model bath until the armihila- 
tion processes drop below the Hubble expansion rate. 
After this point in time the abundance of dark maffer 
is essentially constant throughout the universe. This 
constitutes the weakly interacting massive particle 
(WIMP) paradigm. 

In a recent paper ||T1 the authors revived an alter¬ 
native mechanism Eia for achieving fhe observed 
dark maffer relic densify. Insfead of using 2^2 
annihilation processes they assume that a dominant 
3^2 number-changing process occurs in the dark 
sector involving strongly interacting massive parti¬ 
cles (SIMPs). This process reduces the number of 
dark particles at the cost of heating up fhe secfor. 
However, fhe presence of hof dark matter is problem- 
afic for sfrucfure formation, which means that at the 
time of freeze-ouf, fhe dark maffer particles must to 
be in thermal equilibrium with the standard model 
ones. This, in turn, requires small couplings between 
the dark and the standard model sectors. In this way 
the energy from fhe dark secfor can be fransferred fo 
the standard model via scattering processes. 

The coupling between the two sectors allows for 
direcf and indirecf defection, while fhe large self- 
inferacfions can play a role in sfrucfure formafion, 
by solving fhe core vs. cusp problem IH. Compared to 
the WIMP paradigm, where the dark matter particles 
typically are believed to be around the TeV scale, this 
model can yield dark matter particles with masses 
around a few 100 Me Vs. This is an inferesfing alfer- 
nafive fo fhe WIMP paradigm given fhe facf fhaf cur- 
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renf experimenfs are puffing subsfanfial consfrainfs 
on fhis old paradigm. These consfrainfs can, how¬ 
ever, be alleviafed or even be complefely offsef wifhin 
fhe recenfly proposed safe dark matter paradigm 13 • 

A follow-up paper ||6| infroduced a realization 
of fhe SIMP mechanism based on an underlying 
sfrongly coupled secfor described via chiral perfur- 
bafion fheory. In fhis sef-up fhe pions consfifufe 
fhe dark maffer particles and a key role is played 
by fhe time-honoured Wess-Zumino-Witten (WZW) 
term 130. The WZW term is non-vanishing in the¬ 
ories where the coset space of fhe symmefry break¬ 
ing paffern has a non-frivial fifth homofopy group. 
This topological ferm infroduces a 5-poinf pion inter¬ 
action, making if an ideal candidate for fhe 3 —> 2 
annihilafion process. In QCD, for example, fhe ferm 
describes fhe annihilation of fwo kaons into three pi¬ 
ons. In this paper we shall only be concerned with the 
symmetry breaking pattern SU{2Nf) —> Sp{2Nf) for 
Nf = 2. The simplesf realizafion of fhis breaking paf¬ 
fern comes from an underlying Sp(2) = SU{2) gauge 
group, buf in general if can be realized for any SpfNc) 
gauge group. The acfual paffern of chiral symmefry 
breaking depends on fhe number of flavours, colors 
and maffer represenfafion. A comprehensive sfudy 
of fhe conformal window for SpfNc) gauge groups 
can be found in I flOl . Laffice simulations have fur- 
fher demonsfrafed fhaf such an underlying d 3 mamics 
fruly leads fo fhe expecfed paffern of chiral symme¬ 
try iHl with the spectrum of the composite spin-one 
resonances first computed in iTTlffT^ . 

The computations performed in |(6| make use of the 
first non-vanishing order in the chiral expansion for 
the 3 —> 2 and 2 —> 2 processes. For the 3 —> 2 annihi¬ 
lation process it is one order higher than the related 
2 —> 2 self-interaction process. We will demonstrate 
that it is important to analyse the physical results via a 
consistent next-to-leading (NLO) and next-to-next-to- 
leading order (NNLO) chiral perturbative treatment. 
We will, in fact, show that the leading order analysis is 
phenomenologically unreliable because it is outside 
fhe range of convergence. The imporfanf higher order 
correcfions subsfantially increase fhe range of conver¬ 
gence of fhe fheory, and fherefore fhe phenomeno¬ 
logical reach. Wifhin our confrollable analysis we 
discover fhaf higher order correcfions increase fhe 
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tension with respect to the phenomenological con¬ 
straints making it hard for the SIMPlest realisation to 
be phenomenologically viable. 


II. CONSISTENT SETUP 


Chiral perturbation theory is an expansion in pow¬ 
ers of fhe pion mass and momenfum. Throughouf 
fhis paper we will use fhe sfandard ferminology LO, 
NLO and NNLO when referring fo 0{p^), 0{p^) and 
0{p^), respecfively Furfhermore, beyond leading or¬ 
der if is well known fhaf a consisfenf chiral perfur- 
bafion analysis requires fhe infroducfion of a number 
of low-energy consfanfs fhaf uniquely identify fhe 
underlying sfrongly coupled fheory. The NLO low- 
energy consfanfs, relevanf for fhe 2 —> 2 process, are 
denofed here by N-LECs. Af NNLO bofh processes 
confain low-energy consfanfs which we denofe by 
NN-LECs. 

Since fhe 3 —> 2 process emerges nafurally af fhe 
four-derivative level, if does nof have a leading order 
confribufion. On fhe ofher hand fhe 2^2 process 
does have leading order confribufions, infroducing a 
pofenfial mismafch in power counfing befween fhe 
two processes. In our sef-up fhe solufion fo fhe Bolfz- 
marm equation (which depends on fhe fhermally av¬ 
eraged cross section for fhe 3^2 process) refurns 
a value for fhe pion decay consfanf, as a function of 
fhe pion mass, such fhaf fhe expecfed relic densify is 
obfained. These values are subsequenfly subsfifufed 
info fhe cross secfion for fhe 2 —> 2 self-inferacfion, 
which is consfrained by observafions of e.g. fhe bullef 
clusfer ITU . Technically, fhe mismafch arises because 
fhe two processes, unless computed to the same or¬ 
der, disagree on the definition of the physical pion 
mass m-n and decay constant /„ which at NLO can be 
written as: 
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Here m and / are the bare parameters from the lead¬ 
ing order Lagrangian, a, and bf are combinations of N- 
LECs and L is a log term defined in equation | [T^ . The 
amplitude should, however, be expressed in terms of 
the physical quantities by rewriting the bare quanti¬ 
ties in terms of the physical ones. Erom this argumen¬ 
tation it is clear that a consistent calculation requires 
both processes to be calculated at the same order. 

In practice it means that one has to go beyond lead¬ 
ing order when computing the 2^2 process, and 
this is what we will do here. The additional benefit of 
extending the computation of the 2^2 processes to 
NLO is that we will also be able to estimate, and ex¬ 
tend, the convergence of the perturbative expansion 
towards the physically relevant regime. To check the 
stability of the perturbative expansion we also esti¬ 
mate the NNLO corrections. Because we are inter¬ 


ested only in a preliminary estimate of the NNLO 
corrections we retain the direct contributions from 
the loop diagrams and neglect the finite contributions 
stemming from the NN-LECs. We will show that the 
NNLO results are already much closer to the NLO, 
than the NLO are to the LO for a much wider range 
of pion masses. 


III. SCATTERING AND RELIC DENSITY 

The Boltzmarm equation for the 3 —» 2 process 
reads 


h + 3Hn = -(n^ - n^net,){av^)z^ 2 , (3) 

where n{t) is the pion number density and H is the 
Hubble constant. As argued in [T| contributions from 
the 2^2 armihilation into standard model particles 
can be neglected, because they are sub-dominant. To 
solve the differential equation we rewrite it in terms 
of the dimensionless quantities Y = n/s and x = m^/T 
with s being the entropy and T the temperature. Eor 
the entropy and matter density we used the defini¬ 
tions 


P = s = (4) 

where ge and hg are the effective degrees of freedom 
contributing to the entropy and density, respectively. 
In this notation, the differential equation becomes 


dx 


47i5y. 


ml 


91125G x5 


■(Y3 - Y^Y,Mav^h^2 


(5) 


where gt{x) combines information from both ge(T) 
and he{T) into a single function for the effective de¬ 
grees of freedom (the function is different from the 
function used in the 2^2 case). The equilibrium 
function Y^^ can be written as 


ye^ = 


ASNnX^ 

An^ge{m„lx) 


K2{x), 


( 6 ) 


where K 2 {x) is the modified Bessel function of the 
second kind and Nn is the number of pions. The 
thermally averaged scattering cross section used in 
reads 


{01?h^2 = 


75 Vs mlNl 
512n7x'^ /lo ' 


(7) 


where Nc is the number of colors and our definition 
of the pion decay constant /„ differs by a factor of 
one-half compared to the original paper. The value 
of fn is chosen in such a way that the solution to the 
Boltzmarm equation gives the expected relic density. 
An estimate of the higher order corrections to the 
thermally-averaged cross section can be found in the 
appendix. 
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FIG. 1: Dashed lines belong to the left axis and solid lines to the right axis. The dashed orange line is the solution winZ/n 
to the Boltzmann equation and the dashed horizontal line is the upper perturbative limit m-nl fn = 4ti. The two solid lines 
are the cross section for the 2 —> 2 self-interactions at LO (blue) and NLO (orange). The band on the orange line is the 
uncertainty from the low-energy constants (N-LECs). The solid grey band is the upper limit on the self-interactions. 


A. Pion-pion Scattering 


To determine the pion-pion scattering amplitude 
we make use of the non-linearly realized chiral La- 
grangian. We adopt the notation of ITSl together with 
the results from the NNLO. In fact, we have inde¬ 
pendently performed the NLO computations for the 
breaking pattern SU{2Nf) —> Sp{2Nf) and shown that 
the earlier published results IT^ must be amended. 
To keep the work self-contained we briefly outline the 
computational setup. 

Let G be the global flavor symmetry of the vector¬ 
like fermions transforming according to a given rep¬ 
resentation of the underlying gauge d 3 mamics, and let 
H be the stability group after spontaneous symmetry 
breaking. The associated Goldstone boson manifold 
G/H is parametrized by 


u = expf— 

lV2/ 


( 8 ) 


where X“ are the broken generators, normalized as 
(X'’X*’) = 6®^ where (•) denotes trace in flavor space. 
The quantity u transforms as 


u —> guh^ = hug^, 


(9) 


with g £ G and ti e H. Here ^ is a space-time indepen¬ 
dent element of G while h is a space-time dependent 
element of H that must satisfy the constraint equation 
ID- The chiral Lagrangian quantities that transform 
homogeneously under the stability group H are 

Wf, = - irfi)u - M(d,, - (10) 

= u'xu' + ux'u. (11) 

In the first expression and 1^, are external currents 
that can be set to zero in the computation of the pion- 
pion scattering amplitude. In the second expression 


y is a spurion field that formally ensures chiral sym¬ 
metry invariance at every step in the computation. 
Only at the very end it is replaced by its expectation 
value X - which explicitly breaks chiral symme¬ 
try. The pion mass is then written as 
where Bq is related to the underlying chiral conden¬ 
sate and niq is the new quark mass. In this notation 
the LO Lagrangian is 

£2 = ^(u^u^‘+ x+), (12) 

and the relevant part of the NLO Lagrangian is 


+ L2{UnUv){id‘u") + Lj,{UuU^Uvu'’) 

(13) 

-I- L4,(u^‘u^Xx+} + L5(Mf'w,,X+> 

+ U(x+f + UiX-f + lUixl + X-}- 

Because the NLO Lagrangian represents the most 
general Lagrangian it is possible to absorb the one- 
loop divergences by an appropriate renormalization 
of the low energy constants. We use the modified MS 
scheme where 

U = u.- -^R, (14) 

' 327x2 ^ ' 

with 

X = ^-t-log(47x) - Pe + 1- (15) 

Heree = 4-dandyE = -T'(l)istheEuler-Mascheroni 
constant. It should be noted that the renormalized 
coefficients U will depend on the energy scale p in¬ 
troduced by dimensional regularization. To simplify 
the results later on we use the following short-hand 
notation for the log terms arising from the loop dia¬ 
grams. 


L = 7X16 log 



167X2 • 


(16) 
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FIG. 2: Dashed lines belong to the left axis and solid lines to the right axis. The red dashed line is the NNLO solution winZ/n 
to the Boltzmann equation, the orange dashed is the NLO, and the dashed (grey) horizontal line is the upper perturbative 
limit m„//n = in. The three solid lines are the cross section for the 2 —> 2 self-interactions at LO (blue), NLO (orange) and 
NNLO (red). The purple band is the uncertainty from the low-energy constants (N-LECs). The solid grey band is the upper 
limit on the self-interactions. 


At NLO there are three new diagrams contributing 
to the amplitude for pion-pion scattering (see IflSl for 
defails and diagrams) and af NNLO fhere are twelve 
addifional diagrams. The scattering amplifude can 
be written as 

T{s,t,u) = 

+ e^‘’‘^B{u,s, t) + 6‘'^6‘^‘^C(s, f, u) (17) 
-t 5‘"=5'’‘^C{t,u,s) + 6‘"^6'"^C(M,s,f), 


The gauged version of fhis term and its generaliza¬ 
tion to include left and right vector sources can be 
found in IT^ . 

Af NLO fhere is only one diagram for fhe 3 —» 2 
process while af fhe NNLO fhere are fhree more dia¬ 
grams. However, only fwo of fhese diagrams are non¬ 
zero in fhe approximation of vanishing NN-LECs. 

IV. RESULTS 


where is a group fheorefical facfor fhaf depends 
on fhe breaking paffem and fhe number of flavors. 
For fhe case 511(4) —> Sp(4) all relevanf quantifies can 
be found in fhe appendix. 


B. Topological term 


The topological term can be written in a compact 
local form as a five-dimensional infegral I3-l9l known 
as fhe Wess-Zumino-Wiffen (WZW) ferm 


Nc 


’WZW 


2407t2 


J" da J" d^x e"^‘^*(waWfcWcWdMe), (18) 


were a is fhe fifth spacetime coordinafe. We furfher- 
more redefine fhe fwo quantifies 


( 


ta 


■X“ 


w = exp —— 

lV2/ 

Ua = i{u^daU - udaU"), 


(19) 

( 20 ) 


with da being a five-dimensional derivative. Ordi¬ 
nary Minkowski space is now defined as fhe surface 
of fhe five-dimensional space where a = 1. Furfher- 
more fhe pre-factor of fhe WZW ferm confains direcf 
informafion abouf fhe underlying gauge d 3 mamics. 


We are now ready fo discuss fhe physical resulfs 
bofh af fhe NLO and NNLO. 

Af fhe NLO we refain fhe N-LECs fhaf uniquely 
specify fhe underlying sfrongly coupled fheory. 
These consfanfs are expected fo be a number of or¬ 
der u nify suppressed by a facfor Tiig, defined in equa- 
fion < [T^ . To determine fhe LECs one can eifher use 
experimenfs, when available (as for QCD), or per¬ 
form laffice simulafions when fhe underlying fheory 
is known lUTHTSlI . For QCD fhe N-LECs are known 
ITtI and mosf of fhem are of fhe order O{10~^). In- 
feresfingly fhis is an order of magnifude smaller fhan 
naively expected. Therefore fo estimate fhe N-LECs 
we randomize fheir values using a Gaussian disfribu- 
fionwifh spread 5x10“^ and zero mean. By averaging 
over a large number of samples we arrive af a NLO 
band fhaf reflecfs fhe size of fhe confribufions from 
fhe N-LECs. 

In figure we plof fhe solufion fo fhe Bolfzmann 
equation (dashed orange line) and fhe cross secfion 
for fhe 2 —> 2 self-inferacfions (solid lines) af LO (blue 
line) and NLO (orange line). The band on fhe or¬ 
ange line is fhe estimated effecf from fhe low-energy 
consfanfs. The horizonfal dashed line is fhe upper 
perfurbafive limif nin!fn = Ttt. The cross secfion for 
fhe self-inferacfions is consfrained from above by fhe 
solid grey line. The confrainf from fhe bullef cluster 
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FIG. 3: The shaded regions indicate where the different orders in perturbation theory can be trusted up to around 20%. The 
blue line is LO, orange line is NLO and red line is NNLO. The dots show for, each order, when the next order corrections 
are of the same magnitude. The estimates are obtained neglecting the LECs. 


IT^ is a Inin ^ 1 cm^/g but simulations of halo struc¬ 
tures IITSlITgl suggest that the limit might be closer to 
olm-n < 0.1 cm^/g. This uncertainty is reflected in the 
grey band. 

The model is phenomenologically relevant when 
the dashed orange line is below the perturbative limit 
and the self-interactions are below the solid grey line. 

However, for Nc ~ 2 these two criteria are never 
met simultaneously when the NLO corrections are 
properly taken into account. The results demonstrate 
that the NLO corrections are crucial to establish the 
phenomenological viability of the model. The right 
panel of figure]^ shows that by substantially increas¬ 
ing the number of colors one may still hope to meet 
the phenomenological criteria although in this case a 
more careful use of the large-Nc expansion is needed. 

Because we assume small couplings to the standard 
model, there is also a lower limit on the dark matter 
mass of around 10 MeV. In fact, lighter masses will 
change the effective number of neutrino species 1201 
which is in tension with the Planck data ||2T|. 

Since, for the region of phenomenological interest, 
the NLO corrections are quite sizeable (especially for 
small values of Nc) we also estimate the NNLO correc¬ 
tions. To consistently include these higher order cor¬ 
rections, we re-determined the solution to the Boltz- 
marm equation using the thermally averaged cross 
section 

, ( 21 ) 

at NNLO. The term comes from the loop cor¬ 
rections whereas the last term bn, would contain the 
NN-LECs if we had included them. As it turns out, 
the solution is relatively insensitive to the value of bu, 
assuming it is of order 0(7116). This means that setting 
the low-energy constants to zero in the amplitude for 
the 3-^2 process is a reasonable approximation. For 
the self-interactions we estimate the effects from just 


the N-LECs in the same way as before. 

In figure 0 we show the solution to the Boltzmann 
equation (dashed red line) and the 2 —> 2 scatter¬ 
ing cross section (solid red line) at NNLO. There are 
two important points to make, namely that (a) the 
solution nin/fn decreases, making the investigation 
more amenable to a perturbative analysis, while (b) 
the cross section seems to increase further, strength¬ 
ening the tension with the self-interaction constraints. 
If the solution did not decrease one would have seen 
a more dramatic increase in the NNLO cross section. 
This underlines the importance of a consistent calcu¬ 
lation. Despite the small difference between the NLO 
and NNLO cross sections, we observe a more promi¬ 
nent effect from the LECs as indicated by the purple 
band. There are two reasons for the increased uncer¬ 
tainty. It is partly due to the fact that the N-LECs enter 
more frequently in the NNLO result, but also because 
the squared norm of the NLO amplitude enters in the 
NNLO cross section. Because this quantity is strictly 
positive, the average of a Gaussian distribution devi¬ 
ates from zero. 

In figurel^we provide regions of validity where the 
different orders in perturbation theory can be trusted 
up to around 20%. The dots show where each order 
breaks down because the corrections from the next 
order are of equal magnitude. For the NNLO case 
the values are rough estimates from extrapolating via 
the previous orders. In the grey region the theoreti¬ 
cal uncertainty of the NNLO expansion is more than 
20%. We therefore estimate that for Nc = 2 pertur¬ 
bation theory at NNLO is reasonable up to around 
30 MeV. However, besides the aforementioned con¬ 
straint coming from the effective number of neutrino 
species for masses below 10 MeV, this region is ex¬ 
cluded by the constraint on the self-interactions. Even 
by allowing the NNLO treatment to be used beyond 
this point, the model cannot easily meet the phe¬ 
nomenological constraints. It is clear from the figure 
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that the claimed phenomenologically relevant region 
161 is far beyond LO perturbation theory. 

For the case of Nc = 16 one has confrol up fo around 
110 MeV af NNLO, buf here fhe fheory is sfill in fen- 
sion with the constraints on self-interactions. How¬ 
ever, after this point the theory meets the constraints 
for plausible combinafions of low-energy consfanfs. 

One can introduce additional breaking of fhe 
flavour symmefry which, in principle, can reduce fhe 
2^2 scaffering amplifude but it will not change 
the range of convergence for the different perturba¬ 
tive orders. 


V. CONCLUSION 


SLr(4) —> Sp(4) pattern of chiral symmetry breaking is 
^abcd ^ ^ ^ad^bc^^ (^ 2 ) 

with {a, b,c,d] = 1,..., 5. The functions B(s, f, u) and 
C(s, t, u) evaluated on the kinematical threshold read: 

LO - All terms should be multiplied by jf^. 

B(s, t,u) - 1 
B{u,s,t) = -1 
C(s, t,u) - 0 
C{u,s,t) = 0 


By consistently using chiral perturbation theory we 
studied the phenomenological viability of an interest¬ 
ing class of composite dark matter models. In these 
models the relic density of the dark pions is achieved 
via 3 —> 2 number-changing processes making use of 
the WZW term. We determined both the 3^2 and 
the 2 —> 2 processes to the NLO and NNLO order in 
the chiral expansion and showed that higher order 
corrections substantially affect the LO result of |6] for 
the 2^2 processes. At the NLO and NNLO we have 
shown that the SIMPlest models ||6l, with a moderate 
number of underlying colors, are at odds with phe¬ 
nomenological constraints. However, for sufficiently 
many colors, or via additional breaking of the flavor 
symmetry, it is still possible to find small regions in 
the parameter space consistent with phenomenologi¬ 
cal constraints. 

Our results are also of immediate interest for differ¬ 
ent realizations of composite dark matter 1221 and/or 
composite (Goldstone) Higgs Il23ti25l featuring the 
same pattern of chiral symmetry breaking. Further¬ 
more first principle lattice simulations of the under¬ 
lying d 5 mamics ITTtiT^ are able to provide crucial 
information on the low-energy constants l26l . 
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Appendix A: Scattering Amplitude 


The cross section used for the plots has been eval¬ 
uated at the scattering threshold where s = 4m^ and 
t = w = 0. In this limit the cross section reads 


02^2 


|TP 

128nN^m^ 


(Al) 


with T the amplitude of the process. The group 
theoretical factor in the scattering amplitude for the 


NLO - All terms should be multiplied by 

B{s,t,u) = \{2ai + 32 a 4 -i- /ttis) 

B{u,s,t) = (ai -I- 4 a 2 + 16 a 3 - 7 ii 6 ) 

C(s, t, u) = liApi + 16/32 + 64 jS 3 + IVnu) 
C{u,s,t) = (^1 + 16/34 + nig) 

NNLO - All terms should be multiplied by 

B{s, t, u) = + (02^16 -I- 16 y 4 + yi 

B(u,s, t) = + coatik, + 64 y 5 + 16^3 + 4 y 2 + Ti 

C(s, t, li) = + (UgTIlg -I- 6465 + 1663 •+■ 462 + b\ 

C{u,s,t) — + (Ugnig -1-1664 -1- 61 


Here we made use of crossing symmetry which dic¬ 
tates that B(t,u,s) = B(s,t,u) anci C(t,u,s) = C(u,s,t). 
The definitions of oji are given by 


0)2 = 16(L;, + U^+ 3Ll + L'g) 

+ 128{L[ +L\-L\+ LD - f L 

C 03 = 

0)4 = 32{Ll + LI-LI+ LD 

- 128(L!i + + LD + f L 

ri - 20 _ 14„2 
CO^ 3 9 ^ 

a;6 - 152(L^ + LJ + L^) - 8L^ 

+ 416(L!i +L'2+LI+ LD - 223l 
"7 = - I 

0)8 = -32(L;, -h + L^) - 2L 


The remaining constants a„ /S,, y, and 6; can be found 
in 1T5]| . For the breaking pattern SU{2Nf) —> Sp{2Nf) 
considered here, we amended the expression for a 2 
given in m where the sign in front of the Tiig term 
should be negative. We furthermore set all LECs at 
NNLO to zero. We choose the renormalization scale 
p2 = 4 jj 22 agree with the energy scale of the process. 
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Appendix B: Thermal Average 


with the following definitions 


The thermally averaged cross section for fhe 3 —> 2 
process af NLO is 






and af NNLO reads 


75 Vs mlN^ 
5127i5x2yio jys' 


(Bl) 


^ , (B2) 


UiD 


5^/5 / 9+ V45 \ 

2887i2 °^\9_ Visj 


(B3) 


7 

967t2 


^ 1.83 • IQ-'^ 


(B4) 


If we had included fhe NN-LECs fhey would appear 
in b^a. 
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